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The thermal properties of the glueballs are studied using SU(3) anisotropic lattice QCD with 
Aat = 6.25, the renormalized anisotropy £ = a s /at = 4 over the lattice of the size 20 3 x N t 
with N t = 24, 26, 28, 30, 33, 34, 35, 36, 37, 38, 40, 43, 45, 50, 72 at the quenched level. To construct 
a suitable operator for the lowest-state glueball on the lattice, we adopt the smearing method, 
providing an illustration of its physical meaning in terms of the operator size. First, we construct 
the temporal correlators G(t) for the lowest ++ and 2 ++ glueballs, using more than 5,500 gauge 
configurations at each temperature T. We then perform the pole-mass measurement of the thermal 
glueballs from G(t). For the lowest ++ glueball, we observe a significant pole-mass reduction of 
about 300 MeV in the vicinity of T c or ma{T ~ T c ) ~ 0.8ma(T ~ 0), while its size remains almost 
unchanged as p(T) ~ 0.4 fm. Finally, for completeness, as an attempt to take into account the effect 
of thermal width T(T) at finite temperature, we perform a more general new analysis of G(t) based 
on its spectral representation. As an ansatz to the spectral function p{uj), we adopt the Breit-Wigner 
form, and perform the best-fit analysis of the temporal correlator as a straightforward extension 
to the standard pole-mass analysis. The result indicates a significant broadening of the peak as 
r(T) ~ 300 MeV as well as rather modest reduction of the peak center of about 100 MeV near T c 
for the lowest ++ glueball. The temporal correlators of the color-singlet modes corresponding to 
these glueballs above T c are also investigated. 

PACS numbers: 12.38.Gc, 12.39.Mk, 12.38.Mh, 11. 15. Ha 



INTRODUCTION 



QCD at finite temperature is one of the most interest- 
ing subjects in the quark hadron physics |IJ g, |J. At 
low temperature, quarks and gluons are confined in the 
color-singlet objects, hadrons, due to the nonperturba- 
tive effect of the strong interactions among them. With 
the increasing temperature, these interactions diminish 
due to the asymptotic freedom, which is expected to lead 
to liberation of quarks and gluons above a critical tem- 
perature T c , forming the new phase, i.e., the quark-gluon 
plasma (QGP) phase. 

At present, we are in front of the experimental possi- 
bility to create the QGP phase in the RHIC project at 
Brookhaven National Laboratory. Hence, much progress 
is desired in the theoretical understanding of QCD at fi- 
nite temperature. To study QCD at finite temperature 
and to understand the QCD phase transition, the effec- 
tive models and the lattice QCD Monte Carlo simulation 
provide useful and complementary approaches. While 
the effective models offer their own insights into the phys- 
ical phenomena and often provide analytical methods to 
estimate related physical quantities, the lattice QCD sim- 
ulation provides a model-independent method of calcu- 
lating physical quantities directly based on QCD. 



The SU (3) lattice QCD at the quenched level indicates 
the deconfinement phase transition of the weak first order 
at the critical temperature T c ~ 260 MeV B , and simula- 
tions with dynamical quarks show the chiral phase tran- 
sition at T c = 173(8) MeV for N f = 2 and T c = 154(8) 
MeV for N f = 3 in the chiral limit J§. 

Above T c , most of nonperturbative properties such 
as color confinement and spontaneous chiral-symmetry 
breaking disappear. Consequently, quarks and gluons 
are liberated, and the tremendous changes are expected 
in the mass spectrum. 

Even below T c , one is interested in the hadronic mass 
shift (|, [| |^] , because it can be one of the most important 
pre-critical phenomena of the QCD phase transition at 
finite temperature (and also at finite density). In this 
respect, from the experimental side for instance, CERES 
Collaboration Q has proposed the high-energy heavy-ion 
collision data, which may indicate the mass shift of p- 
meson, and extensive theoretical efforts have been made 
aiming at understanding of its true implications || . 

The theoretical background to believe in the hadronic 
mass shift is as follows. As the temperature increases up 
close to T c , the inter-quark potential is known to change 
significantly (l0[ O. As far as the heavy quarkonia are 
concerned, for instance J/tp, the change of the inter- 
quark potential may be followed by the changes in the 
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structures of hadrons, which may lead to significant mass 
shifts as a consequence 0. One also expects significant 
mass shifts of the light hadrons, for instance, er-meson 
H , which is due to the partial chiral restoration. 

At finite temperature, due to the absence of the 
Lorentz invariance, the concept of the "mass" becomes 
less definitive, and there emerge two distinct concepts of 
the "screening mass" and the "pole-mass" . The screening 
mass governs the correlations along the spatial direction, 
while the pole-mass governs the correlations along the 
temporal direction. The relation between them is anal- 
ogous to the relation between the Debye screening mass 
and the plasma frequency in electro dynamics fll2[ . Both 
of these two masses are interesting quantities, which can 
reflect the important nonperturbative features of QCD. 
However, we are more interested here in the pole-mass 
than in the screening mass, because the pole-mass is 
closely related to the QGP creation experiments and cal- 
culations mentioned above. In fact, the pole-mass is ex- 
pected to be directly observed as the physical mass of the 
thermal hadrons. 

In contrast to the effective model approaches, which 
have been extensively used to study the thermal proper- 
ties of hadrons |2|, |) , only a few lattice QCD studies have 
been performed for the pole-mass of thermal hadrons so 
far g§ 0, g |l(J. Instead, the screening mass has been 
extensively studied with the lattice QCD Jl|, [Tj], [TJ, [TJ 
by analyzing the spatial correlations of the hadronic cor- 
relators. The reason behind this is the technical difficulty 
of measuring hadronic two-point correlators in the tem- 
poral direction at finite temperature on the lattice. In 
fact, since the temporal extension of the lattice shrinks 
as 1/T at high temperature, the pole-mass measurements 
have to be performed in the limited distance shorter than 
1/(2T), which corresponds to Np/2 = 2-4 near T c in the 
ordinary isotropic lattice QCD fll2fl . 

The severe limitation on the measurement of the tem- 
poral correlation can be avoided with an anisotropic lat- 
tice where the temporal lattice spacing at is smaller than 
the spatial one a s p0| . On the anisotropic lattice, by 
taking a small a t , it is possible to use efficiently a large 
number of the temporal lattice points even in the vicinity 
of T c , while the physical temporal extension 1/T — N t at 
is kept fixed. In this way, the number of available tempo- 
ral data are largely increased on the anisotropic lattice, 
and accurate pole-mass measurements from the temporal 
correlation become possible. 

In this paper, we study the glueball properties at fi- 
nite temperature from the temporal correlation in SU(3) 
anisotropic lattice QCD. The glueballs are special and 
interesting hadrons. They belong to a different class 
of hadrons. Unlike mesons and baryons, which mainly 
consist of valence quarks and anti-quarks, the glueballs 
mainly consist of gluons in a color-singlet combination. 
Their existences are theoretically predicted in QCD as a 
consequence of the self-interactions among gluons, and 
numbers of lattice QCD calculations have been per- 
formed aiming at investigating the properties of glueballs 



p2 . 23 1 . In the real world, the glueballs are expected 
to mix with the qq mesons due to the presence of the 
light dynamical quarks. This mixing problem raises a 
difficulty in distinguishing them from ordinary mesons in 
full lattice QCD as well as in experiments. Still, enor- 
mous experimental efforts have been and are being made 
in the pursuit of the glueballs in the real world |24j] . At 
present, /o(1500) and /o(1710) are taken as serious can- 
didates for the 0++ glueball, and £(2230) for the 2++ 
glueball, and these glueball candidates are considered to 
contain substantial fractions of glueball component. In 
this paper, to avoid this difficulty of identifying the glue- 
ball, we adopt quenched lattice QCD as a necessary first 
step before attempting to include the effect of dynami- 
cal quarks in future. It is worth mentioning here that, 
even without dynamical quarks, quenched lattice QCD 
reproduces well various masses of hadrons, mesons and 
baryons, as well as the important nonperturbative quan- 
tities such as the string tension and the chiral condensate. 

In quenched QCD, the elementary excitations are only 
glueballs in the confinement phase below the critical tem- 
perature T c ~ 260 MeV. From the lattice QCD stud- 
ies at quenched level |2l], [23|| , the lightest phys- 
ical excitation at zero temperature is known to be 
the scalar glueball with J PC = ++ and the mass 
m G (0 ++ ) ~ 1500-1700 MeV, and the next lightest one 
is the tensor glueball with J PC — 2 ++ and mc(2 ++ ) ~ 
2000-2400 MeV. These light glueballs are expected to 
play the important role in the thermodynamic properties 
of quenched QCD below T c . Our aim is to understand 
the thermodynamic properties of quenched QCD below 
T c from the view point of the thermal properties of the 
lowest-lying ++ and the 2 ++ glueballs. 

In our previous paper, only the numerical result of the 
pole-mass reduction of the ++ glueball at finite temper- 
ature was presented ]lq| . In the present paper, we go into 
more detail. We reformulate the pole-mass measurement 
itself from the viewpoint of the spectral representation of 
the temporal correlator. In addition to the ++ glueball, 
we also consider the 2 ++ glueball including the decon- 
finement phase above the critical temperature T c as well. 
Furthermore, based on the spectral representation, we 
attempt to take into account the effects of the possible 
appearance of the thermal width of the bound-state peak 
at finite temperature by proposing the best-fit analysis 
of the Breit-Wigner type. 

The contents are organized as follows. In Sect. ||, we 
begin with a brief review of the spectral representation 
of the temporal glueball correlator. We then introduce 
the smearing method as a method to enhance the ground- 
state contribution in the glueball correlator, and consider 
the physical meaning of the smearing method. We also 
give a prescription how to make a rough estimate of the 
glueball size based on the informations from the temporal 
glueball correlators with the smearing method. Sect. Ill 
is devoted to the brief descriptions of the lattice QCD 
action, its parameters, the determination of the scale, 
and the critical temperature T c in our anisotropic lat- 
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tice. In Sect. [V, we construct the temporal glueball 



correlators based on SU(3) lattice QCD, and perform 
the pole-mass measurements assuming that the thermal 
width of the peak in the spectral function is sufficiently 
narrow. We study the low-lying ++ and 2 ++ glueballs 
both below and above T c . This procedure itself has been 
widely used in the standard mass measurements of vari- 
ous hadrons in the lattice QCD Monte Carlo calculations 
at zero temperature. It was also adopted by [l4| |l6| in 
the pole-mass measurements of various mesons at finite 
temperature. The section is closed with the comparison 
of our results with the related lattice QCD results on the 
screening masses and the pole masses at finite temper- 
ature. In Sect. |y|, for completeness, we perform a more 
general new analysis of the temporal correlators of glue- 
balls as an attempt to take into account the effects of 
the possible appearance of the non-zero thermal width of 
the bound-state peak at finite temperature. We adopt 
Breit-Wigner ansatz for the form of the spectral func- 
tion to extract the center and the thermal width of the 
ground-state peak through the spectral representation. 
This procedure itself is a straightforward extension to 
the procedure adopted in the pole-mass measurements, 
which seems to be widely applicable to var ious temporal 
correlators of thermal hadrons. Sect. VI is devoted to 
the summary and concluding remarks. 



II. THE SMEARING METHOD IN QUENCHED 
SU(3) LATTICE QCD 

A. The spectral representation 

Using the glueball operator 4> (t,x), we consider the 
temporal correlator of the glueball as 



G(t) = u(t)m), 



(i) 



where <fi(t) is the zero-momentum projected operator de- 
fined as 

<P(t) = 0o (t) - fao), 0o(*) = -jJ-X>(*.2). ( 2 ) 



with N s = N x N y N z being the number of sites in the 
spatial submanifold. The angular bracket (•) denotes the 
statistical average. The subtraction of the vacuum ex- 
pectation value is necessary for the ++ glueball. The 
summation over x is performed for the zero-momentum 
projection, i.e., the total momentum is set to be zero. 
The glueball operator cj>o(t,x) has to be properly con- 
structed so as to have a definite spin and parity in the 
continuum limit. For instance, the ++ and 2 ++ glueball 
operator are defined as J2f], H^, |25|, ||(| 



4> {t,x) = ReTx(Pi 2 (t,x) + P 2 3(t,£) + P3i(t,x)) ,(3) 
cf>o (t, x) = ReTr (2P 12 (t, x) - P 23 (t, x) - P 31 (t, jg)) , 



respectively, where Pij(t,x) denotes the plaquette oper- 
ator in the i-j plane. 

To consider the spectral representation, we express the 
temporal correlator G(r) in the canonical operator rep- 
resentation as 



G(r)=^)- 1 tr(e-^(r)0(O) 



(4) 



where (3 — l/T denotes the inverse temperature, Z(f3) = 
tr(e~^ H ) the partition function, and H the Hamilto- 
nian of QCD. The field operator <j>{r) is represented 
in the imaginary-time Heisenberg picture as 4>(t) = 



<H0)e 



-tH 



By using the identity G(r) = G{0 — r) 



for the manifest periodicity, we obtain the spectral rep- 
resentation as 



{n \<p\ m) 
~2ZW 



■ exp 



-a- 



,E„ 



En 



(5) 



x cosh 



dui 
2 



{T-P/2)(E n -E m ) 
7r2sinh /W2 V VM/ ') 



where E n denotes the energy of the n-th excited state |n) . 
In this paper, |0) denotes the vacuum, and |1) denotes 
the lowest glueball state. Here, p(oj) denotes the spectral 
function defined as 



-0E„ 



(6) 



x2tt (s {uj - E n + E m ) -6(co-E m + . 

Note that p(to) is positive for positive to, and negative for 
negative ui. Due to the bosonic nature of the glueballs, 
p(u>) is odd in u>. 

The spectral function p(ui) is the residue of the glueball 
correlator in the Fourier representation as 

Jo J-oo 2n m - w 

where lji = —— denotes the Matsubara frequency for 

bosons. It is also the residue of the retarded and ad- 
vanced Green functions Gr(ui) and Ga(w) as Q 

du>' p(w') 



Gr/aH 



2n uj — lu' ± it 



(8) 



(+ for the retarded and — for the advanced Green func- 
tions.) Hence, the information of the physical observables 
such as the mass and the width can be obtained by pa- 
rameterizing the spectral function p(u>) properly and by 
performing the best-fit analysis to the lattice QCD Monte 
Carlo data of the temporal correlator G{t) through the 
spectral representation Eq. (||). 
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B. The smearing method 

For the study of the low- lying glueballs, it is impor- 
tant to extract the contributions of the low-lying states 
in lattice QCD Monte Carlo simulations. The smearing 
method is one of the most popular numerical techniques 
to enhance the contributions from the low- lying states. 
Here, we give a brief review of the smearing method in 
SU(3) lattice QCD in the zero-temperature limit. In this 
limit (3 = 1/T — > oo, the spectral function p(uj) and the 
correlator G(t) are simply expressed as 

p(u) = ^A n {s{u -E n )- 8{u + £„)), (9) 



G(t) 
G(0) ~ 

with A n = | (r 



\0)\ /Z and 

An 



\( n \m\ 2 



(10) 



Note that C„ is a non-negative real number with < 
C n < 1 with the normalization as C n — 1 . It will be 

referred to as the overlap between the state cf>\0} and the 
n-th asymptotic state |n). 

In general, to measure the ground-state mass from 
G(t), one seeks for the region where contributions from 
excited states almost die out, leaving only the ground- 
state component as G(t)/G(0) ~ G\ exp(— Ext). In this 
region, one may perform the best-fit analysis with 



C exp(— mat) 



(11) 



to obtain the ground-state mass tuq — E\ and the 
ground-state overlap G = C\. (Here, as was mentioned 
in the previous subsection, |1) denotes the state of the 
lowest-lying glueball, while |0) denotes the vacuum state. 
Because </>(i) is the zero-momentum projected operator, 
defined in Eq. (|2|), E n denotes the mass of the n-th state.) 
In principle, such a region always exists for enough large 
t in the zero-temperature case. However, in practice, 
it is difficult to use such large t in lattice QCD Monte 
Carlo calculations, since the correlator G(t) decreases ex- 
ponentially with t and becomes so small for the large 
t that it is comparable to its statistical error. Hence, 
to measure the mass of the ground-state, it is essential 
that the ground-state overlap C\ should be sufficiently 
large. In the case of the glueballs in quenched SU(3) lat- 
tice QCD, the ground-state overlaps C\ of the operators 
given in Eq. (|^) are quite small, as long as they are con- 
structed from the simple plaquette operators Pij(t,x). 
In this case, due to the considerable contributions from 
excited states, the measured mass with Eq. ( |ll|) always 
behaves as if it were much heavier. The difficulty of this 
small overlap originates from the fact that the plaquette 
operator has a smaller "size" of 0(a) than the physical 



size of the glueball pq| . This problem becomes severer 
near the continuum limit. Hence, for the accurate mass 
measurement, it is necessary to enhance the ground-state 
contribution by improving the glueball operator. This is 
achieved by making the operator have approximately the 
same size as the physical size of the glueball. The smear- 
ing method provides an iterative procedure to generate 
such an extended operators, which is referred to as the 
smeared operators ^8|, p9[ . 

The smeared operator is obtained by replacing the 
original spatial link variables Ui (s) in the plaquette oper- 
ators by the associated fat link variables u[ n \s). Start- 
ing from uf*\s) = Ui(s), the (n + l)-th fat link variable 



U^ n+1 \s) e SU(3) is defined iteratively from U^' l '(s) so 
as to maximize 



An), 



ReTr 



(V(«+i) (s) y/»)t (s) ) 



(12) 



with 



V} n \s) 



(13) 



where U^(s) 



aU^(s) 

£ U^(s)ut ) (s±j)U^(s + i), 

= Ul n)i (s - ft). Here, a G K is referred 
to as the smearing parameter, which controls the speed 
of the smearing. The summation index j runs over only 
the spatial directions to avoid the artificial nonlocality 
in time. A schematic illustration of V i (s) is shown in 
Fig. 0. U^ n+1 \s) is the closest SU(3) element to v} n \s). 
It is easy to see that Vf^ 1 (s) and Jj\ n+X ' (s) hold the same 
gauge transformation properties with U^ n \s). Hence, 
the smearing method respects the gauge covariance. As 
the physically extended glueball operator, we adopt the 
n-th smeared operator, i.e., the plaquette operator which 
is constructed with the n-th fat link variables u\ n \s). 



A 



A 
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FIG. 1: The schematic illustration of V} n) (s). The thick line 



represents the variable V i (s). Thin lines represent the link 
variables u\ n \s). 



C. The physical meaning of the smearing method 

Here, we estimate the size of the n-th smeared oper- 
ator and consider the physical meaning of the smearing 
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method. To obtain the direct correspondence with the 
continuum theory, we adopt the lattice Coulomb gauge, 
which is defined so as to maximize 



(14) 



s fe=l 



with respect to SU(3) gauge transformations. In the con- 
tinuum limit, we recover the continuum Coulomb gauge 
condition as 



3 

E 

*:=i 



d k A k (x) = 0. 



(15) 



In this gauge, the spatial link- variable Uk(s) becomes 
maximally continuous around 1, and can be safely ex- 
panded in a as U k (s) = e ioj4 *W = 1 + iaA k (s) + 0(a 2 ) 
near the continuum limit. We consider the smearing with 
the smearing parameter a. For convenience, we define 



1 



(16) 



which satisfy p + Aq = 1. As shown in Fig. the 
(n + l)-th smeared spatial gluon field A k (s;n + 1) is 
iteratively described with the n-th smeared gluon field 
A k (s; n) as 



A x (s;n + 1) (17) 
= pA x (s; n) + q (^A v (s;n) + A x {s + y;n) - A y (s + x;ri) - A y (s — y;n) + A x (s -y;ri) + A y (s + x - y;n) 

+A z (s; n) + A x (s + z; n) — A z (s + x; n) — A z (s — z; n) + A x (s — z; n) + A z (s + x — z; n) ^ 

I 



in the leading order in a. The expressions for A y (s; n+ 1) 
and A z (s; n + 1) are similar. After some rearrangements, 
we obtain in the leading order in a 



d n A x (s;n) 
= A x (s;n + 1) - A x (s;n) 

= q^dyA x (s;n)-d y A x (s-y;n) 
+d z A x (s] n) - d z A x (s - z\ n) 
-d x A y (s;n) + d x A y {s - y;n) 
-d x A z (s;n) + d x A z (s - z;n) } 

= q { (d*d x + d^d y + dfd,) A x (s; n) 

3 

-c^dfA fc ( S ;n)}. 



(18) 



k=l 



Here, d k and denote the forward and backward dif- 
ference operators along fc-axis (k — x,y,z), respectively, 
which are defined as 



d k f( S ) EE f( S + k)~ f( S ), d£f(x) EE f( S ) - f( S - k), 

(19) 

respectively. Here, k denotes the unit vector along the 
fc-axis. By applying the continuum approximation both 
with respect to s and n, the difference operators d n , d k 
and d k are then replaced by the derivatives a n d/dn, 
a s d/dx k and a 8 d/dx k , respectively. Here, n ee na n is the 
continuum parameter corresponding to n, and a n 6 R is 
formally introduced as a small constant interval in the n- 
direction, although the final result does not depend on a n 
at all. In the continuum Coulomb gauge using Eq. (15), 



Eq. ( |l8| ) reduces to the diffusion equation 
d 

a n —Ai(x;n) = DAAi(x-n) for i = x,y,z, (20) 
on 



where D — 



works as the diffusion parameter. 



By solving this equation, the n-th smeared gauge field 
Ai (x, n) can be expressed by a linear combination of the 
original gluon field Ai(x) as 



Ai(x;n) = d 6 yK{x-y;n)Ai{y) (21) 



K(x;n) = j exp . 



where the Gaussian extension p is defined as 



2VDn = 2ac 



(22) 



In this way, the n-th smeared gluon field physically corre- 
sponds to the Gaussian-distributed operator of the orig- 
inal gluon field. After the linearization of the gluon 
field, the n-th smeared plaquette is also expressed as the 
Gaussian-distributed operator. Here, the Gaussian ex- 
tension p can be regarded as a characteristic size of the 
n-th smeared operator. 

The expression of p in Eq. (^2|) explains the physical 
roles of the two parameters, a and n. The smearing 
parameter a plays the role of controlling the speed of 
the smearing, while, for each fixed a, n plays the role of 
extending the size of the smeared operator. Note that, 
for the larger a, the speed of smearing is the slower. 

The smearing method was originally introduced to 
carry out the accurate mass measurement by maximizing 
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the ground-state overlap j28| . However, we will use it to 
give a rough estimate of the physical glueball size. In 
fact, once we find the pair of the smearing parameters n 
and a, which achieves the maximum ground-state over- 
lap, the physical size of the glueball is roughly estimated 
with Eq. ©. 

With Eq. (^l|), we try to express the n-th smeared 
++ glueball operator in terms of the original gluon field 
A a ^{x). Near the continuum limit, the zero-momentum 
projected ++ glueball operator $gb(0 ++ ;«) consisting 



of the smeared link- variable (s) is expressed as 

$ GB (0+ + ;n) = i J d 3 xG a lJ {x,t-n)G a l] (x,t-n), (23) 

where G^ (x, t; n) is the field strength associated with the 
smeared gluon field Af (x, t;n). Here, t can be regarded as 
a fixed parameter in this argument, and will be omitted 
hereafter. By inserting Eq. < I ) gb(0 ++ ) is expressed 
with the original gauge field Af(x) as 



*gb(0 ++ ; ?1 ) = i J Sx J d 3 yd 3 z (J^K{x - y;n)A 3 (y) - ^-K{i 



x-y;n)Ai(y) 



d d 
—K(x- z;n)Aj(z) - -^—K(x- z\n)Ai(z) 



(24) 



d 3 yd 3 z 
\ J ir 3 / 2 p 3 
d 3 yd 3 z 



2 
V 
2 

V J TT 3 / 2 P 5 



exp 



dxi 
2p 2 



dXj 

—Ai(y)—Ai(z) + 
dyj dz 3 



3 - yy „ ; | exp 



(y 



2p 2 



Ai{y)Ai{z) 



Here, "• • • " represents the terms which are more than 
quadratic in the gluon field, and are dropped off. In this 
derivation, partial integrations and the Coulomb gauge 
condition have been used. We expect that <&gb (0 ++ , n) is 
associated with the best pair of the smearing parameters 
n and a and works as an approximate "creation" operator 
of the ++ glueball, which creates an approximate single 
glueball state from the vacuum. 

D. The relation to Bethe-Salpeter amplitude 

In the previous subsection, we proposed one possible 
method to give a rough estimate of the glueball size. In 
this subsection, we consider the relation between our es- 
timate of the glueball size and the other estimate from 
the Bethe-Salpeter (BS) amplitude, which was adopted 
inRefs.Q. 

The size of the glueball is actually a nontrivial quan- 
tity. Although the charge radius of the glueball can be 
formally defined, its electric charges are carried by the 



quarks and anti-quarks, which play only of secondary 
roles in describing the glueball state in the idealized limit, 
since the glueball by its nature does not contain any va- 
lence contents of (anti-)quarks. 

In order to investigate the size of the glueball, the BS 
amplitude provides a convenient tool. The BS amplitude 
fy(x,y) is expressed as 

* BS (y^) = (0\My)A t (z)\G(P = 0)>, (25) 

which can be obtained as a solution of the BS equation. 
In the non-relativistic limit, it is expected to reduce to 
the glueball wave function in the first-quantized picture. 
This property has been exploited to estimate the size of 
the glueball in Ref.||. 

The connection between these two estimates is pro- 
vided with Eq. (]24|) in the following way. We consider 
< I'gb(0 ++ ; n) associated with the best smearing with the 
Gaussian extension p. Since the ground-state overlap is 
maximized, the ++ -glueball state can be approximated 
as 



|G(P = 0)) ~ $ GB (0 ++ ;n)|0) (26) 



V J 7T 3 /2 p 5 ^ p 2 J v y 2 p 
2 



7r 3/2p5 j y pi J \ 2p 



tfV 3 exp \M7)M0)) 



where we exploit the translational invariance and partial integrations. Similar to Eq. (£4j), the states in Eq. (|2q ) 
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FIG. 2: The inter-quark potential plotted against the inter- 
quark separation r. The solid circles denote the on-axis data 
of the inter-quark potential. The solid line represents the 
result of the best-fit analysis with Eq. (KOh . 



24, 26, 28, 30, 33, 34, 35, 36, 37, 38, 40, 43, 45, 50, 72, which 
correspond to the various temperatures listed in Table |. 
For all temperatures, we pick up gauge field configura- 
tions every 100 sweeps for measurements, after skipping 
more than 20,000 sweeps for the thermalization. The 
numbers N con { of gauge configurations used in our calcu- 
lations are summarized in Table 

We divide the data into bins of the size 100 to reduce 
the possibly strong auto-correlations. However, the auto- 
correlations are not observed so strong in this simulation 
with the 100 sweep separation, since the separation by 
100 sweeps seems already rather large. Unless other- 
wise stated, the statistical errors are understood to be es- 
timated with the jackknife analysis by regarding each bin 
as an independent data point throughout in this paper. 

As for the smearing parameters, which play the impor- 
tant role in extracting the ground-state contribution, we 
have checked that one of the best sets is provided as 



such as \G(P = 0)) and \Ai(y)Ai(z)) are obtained with 
the gauge configuration generated in lattice QCD. Hence, 
this relation would provide a kind of nonperturbative 
construction of the glueball state in terms of the gluon 
field. 

In the constituent gluon picture [32| , the BS ampli- 
tude associated with Eq. (p5l) reduces to the following 
expression as 

^s(f,0).-^(3-J)ex P (-^), (27) 

in the heavy effective gluon mass limit. 

III. ANISOTROPIC LATTICE QCD 
A. Lattice Parameter Set 

We use the SU(3) plaquette action on an anisotropic 
lattice 

S o=^ — E ReTrfl-P^)) (28) 

+§^£ ReTr (i- p * 4(s) )' 

° s,i<3 

where P^ v {s) & SU(3) denotes the plaquette operator 
in the /i-i/-plane. The lattice parameter and the bare 
anisotropy parameter are fixed as /?i a t = 2N c /g 2 = 6.25 
and 7g = 3.2552, respectively, so as to reproduce the 
renormalized anisotropy £ = a s /at = 4 (2Cj. We adopt 
the pseudo-heat-bath algorithm for the update of the 
gauge field configurations in the Monte Carlo calculation. 

To measure the glueball correlators, we per- 
form the numerical calculations on the lattice of 
the sizes 20 3 x N t with various N t as N t = 



a = 2.1, AUr=40, (29) 

in the present lattice calculation for the lowest ++ and 
2 ++ glueballs. As we shall see later, the temporal cor- 
relator G(i)/G(0) and the pole mass m,Q are almost in- 
sensitive to a particular choice of the smearing number 
N smT) as far as N sml ~ 40 with a = 2.1. In fact, in the 
rather wide range as 30 < N Bmi < 50, G(t)/G(0) and mo 
are almost unchanged, and the ground-state overlap C is 
kept almost maximized. (See Fig. |?].) Unless otherwise 
stated, we adopt the smearing with Eq. (j^|) in this paper. 
Only for the accurate measurement, we seek for the best 
value of N smj: with a = 2.1 fixed, although the results are 
almost the same as those with N smx = 40. 

In this paper, we only consider the case with a = 2.1. 
The essential point of the smearing method is to con- 
struct an operator which has the same size with the 
physical glueball size in order to improve the ground- 
state overlap [ p8[ . Hence, a different pair of the smearing 
parameters (a/, -/V s ' mr ) would work as good as the origi- 
nal (a,iV smr ) does, provided that these two pairs of the 
smearing parameters correspond to the same operator 
size. However, we have to take into account the following 
points when picking up the value of a. On the one hand, 
the speed of the smearing is the faster for the smaller 
a, which, however, results in the coarser discretization 
in n in Eq. (^p|). Hence, a has to be enough large so 
as to construct the ground-state glueball operator. We 
have checked that a > 2.0 seems to work well. On the 
other hand, if we adopt the larger value of a, the speed 
of the smearing becomes the slower. Hence, we have to 
scan the wider range of iV smr to maximize the overlap. A 
practical solution is provided by a = 2.1, since it is one 
of the smallest possible values of a, which, in the same 
time, corresponds to enough fine discretization in n in 
Eq.©. 
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FIG. 3: The cubed Polyakov loop Re(P 3 } plotted against 
temperature. The vertical dotted line denotes the estimated 
critical temperature, T c ~ 280 MeV. The error bars are hidden 
inside the symbols. 



B. Determination of lattice spacings 

In order to evaluate the spatial and temporal lattice 
spacings, a s and at, we perform a separate Monte Carlo 
simulation on the lattice of the size 20 3 x 80. Here, we skip 
10,000 sweeps for the thermalization, and pick up gauge 
configurations every 500 sweeps for measurement of the 
static inter-quark potential. The statistical data are di- 
vided into bins of the size 2, where each bin is thought of 
as an independent data point. In Fig. |^, the inter-quark 
potential is shown. We adopt the parameterization of the 
inter-quark potential as 



V(r) =C + err, 



(30) 



where C , A, a denote the offset, the Coulomb coefficient, 
and the string tension, respectively. By fitting the on- 
axis data with V(r), we obtain C = 0.16217(53), A = 
0.06758(42) and a = 0.00883(12). Hence, by assuming 
■\fo = 440 MeV, the lattice spacings are obtained as 

a' 1 = 2.341(16)GeV, a^ 1 = 9.365(66)GeV, (31) 

namely, a s ~ 0.084 fm and at — 0.021 fm. Hence, the 
spatial lattice size of 20 corresponds to 1.68 fm in the 
physical unit. 



C. The critical temperature 

Here, for completeness, we give an estimate of the crit- 
ical temperature T c on our anisotropic lattice, although 
the value of T c ~ 260 MeV with y/a ~ 420 MeV is rather 
established on the isotropic lattice [SJ . 

To estimate the critical temperature T c , we consider 



P = yj d 3 xTr ^Texpi^ dtA Q (x,t)\ (32) 



d 6 xTr 



where V denotes the volume of the space and T the time- 
ordering operation. As is well-known, its thermal expec- 

tation value (P) — — 

Ncont fc=i 



Pk, averaged over lots of 



gauge configurations, is an order parameter of quark con- 
finement associated with the center Z3 symmetry ||^, ^6| , 
where Pk denotes the value of the Polyakov loop P with 
the fc-th configuration. However, to be strict, it is only 
in the thermodynamic limit with the infinite volume that 
Eq. ([52]) can work as the order parameter to estimate T c . 
Note that, since the size of the lattice is limited, tun- 
nelings from one vacuum to another are unavoidable for 
large simulation time even in the Z 3 broken phase. As 
a consequence, due to the associated cancellation of the 
center Z3-phase factor, (P) finally vanishes on the finite 
lattice, which raises a technical difficulty in estimating 
T c from (P). To avoid this, instead of using the total ex- 
pectation value itself, one often analyzes the scattering 
plot, a plot of Pk for each configuration labeled by k. 

Here, for convenience, we analyze the thermal expec- 
tation value of the cubed Polyakov loop as 



conf 



(33) 



k=l 



which provides the equivalent result for the estimate of 
the critical temperature T c with the scattering plot. We 
note that (P 3 ) is not the order parameter of the confine- 
ment phase in a strict sense, since it is invariant under 
the center Z3 transformation. In general, in the symmet- 
ric phase, the expectation value of the order parameter 
vanishes due to the cancellation of the phase factor asso- 
ciated with the symmetry transformations, which never 
happens in the broken symmetry phase. However, in 
the specific case of (P), it is known from the analysis of 
the scattering plot that the norm of Pk for each config- 
uration is already quite small in the symmetric phase, 
due to which its expectation value (P) almost vanishes 
even without involving the cancellation of the center Z3 
phase factor. In contrast, Pk takes a manifestly non- 
vanishing value for each configuration in the broken sym- 
metry phase, which can be used to identify the critical 
temperature T c . 

Such a tendency is quantitatively described with the 
cubed Polyakov loop (P 3 ). In fact, in the plot of (P 3 ) 
versus temperature T, the critical temperature T c is char- 
acterized as the point where it begins to have a non- 
vanishing value. 

In Fig. |[ (P 3 ) is plotted against temperature. By tak- 
ing into account the finite size effects which smoothen the 
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FIG. 4: The susceptibility \ of the Polyakov-loop plotted 
against the temperature. There is a sharp peak around T = 
280 MeV, which indicates the critical temperature T c of the 
phase transition in quenched SU(3) QCD. 



phase transition, we estimate the critical temperature as 
T c ~ 260-280 MeV, which is consistent with the previous 
studies in Ref. @, i.e., T c /^/a = 0.629(3). 

For more accurate determination of the critical temper- 
ature T c , we investigate the Polyakov-loop susceptibility 
X, which is defined as Q 

x = (n 2 ) - (nf (34) 

{RePexp(-27ri/3), argP G [tt/3, it), 
ReP, argP G [-tt/3, tt/3), 

ReP exp (2iti /3) , argP G [-it, -tt/3). 

In Fig. ||, the Polyakov loop susceptibility is plotted 
against temperature. To obtain x, about 3,000 gauge 
configurations are used at each temperature near T c , 
which arc divided into bins of the size 50. The statisti- 
cal errors are estimated with the jackknife analysis. We 
observe that a sharp peak is located between T = 275 
MeV and T = 284 MeV. Hence, we identify the critical 
temperature of our lattice as T c ~ 280 MeV. 



IV. NUMERICAL RESULTS — TEMPORAL 
CORRELATIONS AND POLE-MASSES 

In this section, we construct the temporal glueball cor- 
relator G(t) at finite temperature, and perform the pole- 
mass measurement. We adopt the procedure used in fl3f| 
as the standard one to extract the pole-mass mc(T) of 
the thermal glueball from the temporal correlator G(t) 
at temperature T by using the best-fit analysis with the 
fit function of single-cosh type considered later. Here, we 
regard each glueball as a quasi-particle and assume that 
the thermal width T of the ground-state peak is enough 
narrow in the spectral function p(u>). From the viewpoint 



of the spectral representation Eq. (||) , such a narrowness 
of the peak is necessary to identify the pole-mass as a def- 
inite object and also to justify this standard procedure 
in a strict sense. (In Sect. 0, we will analyze our lattice 
QCD data without assuming the narrowness of the peak 
as a straightforward extension to the current analysis.) 

Suppose that the thermal width is sufficiently narrow, 
i.e., r ~ 0. Then, the spectral function pioS) receives 
the contribution from the ground-state in the following 
manner: 



p(u>) = 2itA (j)(uj — tuq) — S(u> + toq)^ + 



(35) 



where "• • • " represents the contributions from excited- 
states, and A = A(T) and mc = mc(T) are the strength 
and the pole-mass of the glueball peak at temperature 
T, respectively. Note that the appearance of the second 
delta function is due to the fact that the spectral function 
p(u>) is odd in u>, which reflects the bosonic nature of 
the glueball. Corresponding to Eq. (|35| ), by inserting it 
into the spectral representation Eq. (0)7 the ground state 
contribution to the temporal correlator G(t) is expressed 
as 



G(t) 



A 



1 — e -/3m G 



L-tm G + e -{0-t)m a \ + . 



(36) 



Hence, after the ground state contribution is sufficiently 
enhanced, the appropriate fit function for G(t)/G(Q) is 
the following single hyperbolic cosine as 



g{t) ee C (V tmG + e ~ ( ' 3 -* )mG ) 



(37) 



where C and mc are understood as the fit parameters. 
We will refer to C as the ground-state overlap, and to 
the fit function of g(t) as the fit function of single-cosh 
type for simplicity. 

For the accurate measurement, it is required that G(t) 
is dominated by the ground-state contribution in the fit- 
range. (We consider the determination of the fit range 
in the next subsection with an explicit example.) In this 
case, g{t) can properly represent the ground-state con- 
tribution in Eq. (^), and C and ttiq can be extracted 
accurately. In the actual calculations, the complete elim- 
ination of all the remaining contributions of the excited 
states is impossible. Since each hyperbolic cosine com- 
ponent contributes to G(t) positively as is seen in the 
spectral representation Eq. (|), we always have the fol- 
lowing inequalities as 



5(0) < 1, and C 



9(0) 



1 + e~$ mG 



< 



1 



1 



-/3m G ' 



(38) 



Here, g(0) and C take their maxima, if and only if G{t) is 
completely dominated by the ground-state contribution 
in the whole region < t < (3. In this case, we have the 
following equalities instead: 



g(0) = 1, and C 



1 



(39) 
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We thus see that g(0) ~ 1 and C ~ 1/(1 + e~' 3mG ) ~ 1 
can both work as a criterion to determine whether the 
ground-state enhancement is sufficient or not. We will 
refer to g(0) as the normalized ground-state overlap. 
Note that C differs from g(0) only through the factor 
1/(1 + e~P ma ), which is independent of both the smear- 
ing parameters a and N smr . In this section, we examine 
C to determine if the ground-state enhancement is suf- 
ficient or not. We prefer C rather than g(0) because of 
the following reasons, (i) The ground-state overlap C is 
directly obtained from the best-fit analysis with Eq. ( |37| ) 
as one of the best-fit parameters, (ii) In our previous 
work |Q, C has been used for the criterion. However, 
g(0) ~ 1 has the advantage over C ~ 1 in Sect, [v], where 
the counterpart of C cannot be properly defined, whereas 
the counter part of g(Q) can. At any rate, as is seen in 
Tables | and y, they are not so different from each other 
in numerical value, especially below T c . 



(a) 



A. Thermal 0++ glueball 



In Fig. |5 



we show the ++ glueball correlator for 
the low temperature case T — 130 MeV after the suitable 
smearing as Eq. ( p9| ) . We perform the best-fit analysis of 
the single-cosh type Eq. ( |37| ) with the two fit parameters, 
the ground-state overlap C and the pole-mass tuq of the 
lowest thermal glueball. The solid line denotes the result 
of this best-fit analysis. The fit-range is indicated by the 
vertical dotted lines. 



1. The fit range 

We consider the determination of the fit-range. Even 
after the ground-state enhancement is performed, the 
complete elimination of all the contributions of the ex- 
cited states is impossible, especially in the neighborhood 
of t ~ (mod 0). It follows that the equality Eq. @ 
holds only in a limited interval, which does not include 
the neighborhood of t ~ 0. Hence, for the accurate mea- 
surement of the pole-mass toq , we need to find the appro- 
priate fit-range, where contributions from excited states 
almost die out. To this end, we examine the effective- 
mass plot. The effective mass m e s(t) is defined as the 
solution to the following equation as 

cosh ( (t + 1 - N t /2)a t m cff (t)) 
G(t + 1)/G(t) = ^ ' (40) 



cosh I (t - N t /2)a t m c g(t) 



for a given G(t + 1)/G(t) at each fixed t. In Fig. 5 (b), 
we plot t he eff ective mass m e s(t) against t associated 
with Fig. 5 (a). We see that there appears a plateau, a 
region where m e s(t) takes almost a constant value. In 
this region, it is expected that G(t) consists of nearly a 
single spectral component, and, hence, g(t) can properly 
represent the ground-state contribution in Eq. ([36|) . Note 
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FIG. 5: (a) The temporal correlator G(t)/G(0) of the ++ 
glueball at the low temperature (T = 130 MeV). (b) The cor- 
responding cosh- type effective-mass plot. In both figures, the 
solid lines represent the best single hyperbolic-cosine fitting 
in the plateau, which is indicated by the vertical dotted lines. 
The dashed curve in (a) denotes the result of the best-fit anal- 
ysis of Breit-Wigner type performed in the interval indicated 
by the two open triangles. (See Sect. and Eq. ( pp| ) for 
detail.) 



that, if more than a single state contribute with their 
different masses, a nontrivial i-dependence will appear in 
TO G ff(t). In the zero-temperature case, the plateau in the 
effective-mass plot is widely used as a candidate for the fit 
range in the best-fit analysis of single-cosh type Eq. fl37| ) . 
Here, the role of the smearing method is to suppress the 
excited-state contributions as much as possible, so that 
the plateau can be obtained with the smallest possible t. 

In Fig. U (b), the plateau is indicated by the vertical 
dotted lines, which is used to determine the fit-range. 
In this low temperature case T = 130 MeV, the best-fit 
analysis shows C = 0.953(7) and m G = 1482(25) MeV, 
which seems consistent with ttig — 1500-1700 MeV at 
T ~ MeV. We note that the ground-state overlap 
C = 0.95 becomes enough large owing to the suitable 
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FIG. 6: (a) The temporal correlator G(t)/G(0) of the ++ 
glueball at the high temperature (T = 253 MeV). (b) The 
corresponding cosh- type effective-mass plot. The meaning of 
the solid lines in both figures and the dashed curve in (a) are 
the same as in Fig. @. 



smearing. 



In Fig. 6 (a), we show the ++ glueball correlator for 
the high temperature case T = 253 MeV (< T c ) after 



the suitable smearing as Eq. (29). The solid line denotes 
the result of the best-fit analysis of the single-cosh type. 
The fit-range is indicated by the vertical dotted lines, 
which is determined from the plateau in the correspond- 
ing effective-mass plot Fig. 6 (b). We obtain C — 0.89(1) 
and mc = 1289(26) MeV. Owing to the suitable smear- 
ing, contributions from the excited state are seen sup- 
pressed (C ~ 1). 



2. The ground-state overlap 

In Fig. 0, the ground-state overlap C and the pole-mass 
mc for the ++ glueball are plotted against the smearing 
number N smr for the two cases (i) the low temperature 
case T = 130 MeV (circle), (ii) the high temperature case 
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FIG. 7: (a) The ground-state overlap C and (b) the glueball 
mass 77iG plotted against the smearing number N amr for the 
temperatures T = 130, 253 MeV. The solid circle corresponds 
to T = 130 MeV, and the solid triangle to T = 253 MeV. The 
suitable smearing is found to be in a range as 30 < N amt < 
50, where C is kept maximized and mo is kept minimized 
simultaneously. 



T = 253 MeV (triangle) below T c . Here, we fix one of 
the smearing parameters as a = 2.1. We see that, in the 
region N smT < 30, with increasing N smT , the ground state 
overlap C is growing and the pole-mass tug is reducing. 
Then, in the region 30 < iV smr < 50, C is kept max- 
imized and toq is kept minimized, which implies that 
the ground-state contribution in G(t) is maximally en- 
hanced. Beyond this region, i.e., for N sml > 50, since the 
size of the smeared operator exceeds the physical size of 
the glueball, C begins to decline and toq begins to grow 
gradually. 

The maximum overlap and the minimum mass should 
be achieved at almost the same iV smear . In fact, both of 
these two conditions should work as an indication of the 
maximally enhanced ground-state contribution. In prac- 
tical calculations, these two conditions are achieved at 
slightly different N smelLT . However, the numerical results 
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FIG. 8: The ++ glueball mass ma(T) plotted against tem- 
perature T. The solid circle and the cross denote the thermal 
0++ glueball mass obtained with the maximum overlap condi- 
tion and the mass minimum condition, respectively. The ver- 
tical dotted line indicates the critical temperature T c ~ 280 
MeV in quenched QCD. 



on the pole-mass are almost the same. For instance, from 
Fig. [t], the former condition leads to mc(T = 130MeV) = 
1482(25) MeV at 7V smcar = 39 and m G (T = 253MeV) = 
1294(26) MeV at N smcal = 37, while the latter condi- 
tion leads to m G (T = 130MeV) = 1481(25) MeV at 
iVsmear = 41 and m G (T = 253MeV) = 1284(27) MeV 

at Nsmcar = 46. 

As we mentioned before, we divide the statistical data 
into the bins of the size 100 to reduce the possible auto- 
correlations. However, the bin-size dependence is not 
seen so strong even in the vicinity of the critical temper- 
ature T c . Recall that the measurement interval of the 
gauge field configuration is 100 sweeps. Since this inter- 
val is already rather large, the bin-size of 100 seems to 
be sufficient to overcome the auto-correlations. This is 
also the case in the 2 ++ glueball. 



3. The main result on the pole-mass 

From the analysis at various temperatures as T — 130, 
187, 208, 218, 234, 246, 253, 260, 268, 275, 284 MeV, 
we plot the SU(3) lattice QCD result for the pole-mass 
mo (I 1 ) of the lowest ++ glueball against temperature T 
in Fig. |8|. The solid circle and the cross denote the ther- 
mal glueball mass m G (T) obtained with the maximum 
overlap condition and the mass minimum condition, re- 
spectively. We see that they are almost the same. In 
Fig. ||, we observe a significant pole-mass reduction for 
the lowest ++ glueball near T c as 

m G (T) = 1200 ± lOOMeV for 0.9T C < T < T c , (41) 

in comparison with mc{T ~ 0.5T C ) ~ 1480MeV in our 
data or m G (T ~ 0) ~ 1500-1700 MeV in the lattice data 



in pi], |2^, g3|. In this way, we observe nearly 300 MeV 
reduction or about 20 % reduction of the pole-mass of 
the lowest ++ glueball in the vicinity of T c as 

m G (T~0)-m G (T~T c )~ 300MeV, (42) 
m G (T~T c )~0.8m G (T~0). 

It is remarkable that the pole-mass shift obtained here 
is much larger than any other pole-mass shifts which has 
been ever observed in the meson sector in lattice QCD 

PHI- 

We next consider the glueball size. To this end, we seek 
for iV smr , which realizes the maximum overlap condition. 
Such AT smr is listed in Table | for each temperature. By 
using Eq. (^2|), we give rough estimates of the thermal 
glueball sizes p(T) at various temperatures, which are 
shown in Table QL We find 



p ~ 0.4fm, 



(43) 



both at low and high temperature below T c . 

In Table |, we summarize the lowest ++ glueball 
mass m G (T), the ground-state overlap (7 max , correlated 
X 2 /Nuf, the fit-range {ti,t 2 ), the normalized overlap 
g(0), the corresponding smearing number iV smr , the ++ 
glueball size p(T) , and the number of gauge field config- 
urations -/V con f. 



B. Thermal 2++ glueball 



We show, in Figs. |9 (a) and 10 (a) , the 2 ++ glue 



ball correlators G(i)/G(0) for T = 130 and 253 MeV, 
respectively, after the suitable smearings as Eq. (ps|). 
The solid lines denote the results of the best-fit analy- 
sis of single-cosh type. The fit-ranges are indicated by 
the vertical dotted lines, which are determined from the 



in Figs. 9 (b) and 


10 (b) 




From Fig. 3 (a) 


, we find G ~ 0.96 and m G 



2154 



130 MeV, which seems consistent with m G = 2000-2400 
Me V for th e 2 ++ glueball at zero temperature (2l). From 
Fig. " 



10 (a) 



we find G ~ 0.90 and m G ~ 1996 MeV for 
the 2 ++ glueball in the high temperature case T = 253 
MeV (< T c ). 

By combining the results at various temperatures, we 
plot the pole-mass of the thermal 2 ++ glueball against 
temperature in Fig. |ll]. The solid circle and the cross 
denote the pole-masses of the thermal 2 ++ glueball ob- 
tained with the maximum overlap condition and the mass 
minimum condition, respectively. We see again that they 
are qualitatively almost the same. The thermal 2 ++ glue- 
ball shows a tendency to decline by about 100 MeV for 
T < 0.9T C , which is rather modest in contrast to the ++ 
glueball. However, in the very vicinity of T c , it shows a 
sudden reduction of about 500 MeV. 

In Table || we summarize the lowest 2 ++ glueball 
mass m G (T), the ground-state overlap G max , correlated 
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TABLE I: The pole-mass m,G{T) of the lowest ++ glueball at finite temperature T in SU(3) lattice QCD. The temperature T, 
the temporal lattice size Nt, the thermal glueball pole-mass maiT), the maximum ground-state overlap (7 max j the correlated 
X 2 /A r DF, the fit range (ti,ts), the normalized ground-state overlap g(0), the smearing number N B mt, the thermal glueball size 
p(T) with Eq. (^), and the number the gauge configurations iV con f are listed. The best smearing on N sml is determined with 
maximum overlap condition. Above T c ~ 280 MeV, the data such as mo is to be understood as the best-fit parameters. 
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35 
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34 
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0.81(2) 


0.12 
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33 
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1.37 
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0.84(1) 
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( 6,13) 
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6800 
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7700 



% 2 /-/Vdf, the fit-range (t\,t-2), the normalized overlap 
g(0), the corresponding smearing number N sml , the 2 ++ 
glueball size p(T), and the number of gauge configura- 
tions A^ CO nf- 



C. Glueball correlations above T c 

In the deconfinement phase above T c , neither hadrons 
nor glueballs are elementary excitations any more. In- 
stead, quarks and gluons appear as elementary excita- 
tions in the quark gluon plasma. Therefore, we might 
wonder if there is anything interesting in investigating the 
color-singlet modes such as hadrons and glueballs above 
T c . Nevertheless, based on the lattice-QCD Monte Carlo 
studies on the screening mass, it was pointed out that 
some of the strong correlations survive even above T c in 
the scalar and pseudoscalar color-singlet modes such as a 
and pions[^2). Hence, at least in the neighborhood of T Cl 
it may be possible that some of the nonperturbative ef- 
fects survive in the deconfinement phase, which provides 
interesting information of the high-temperature QCD. 

In this subsection, we attempt to examine the tem- 
poral correlator of the glueball-like color-singlet modes 
above T c . We simply apply the pole-mass analysis 
of the temporal correlators on the anisotropic lattice 
with N t = 24,26,28,30,33, which correspond to T = 
390, 360, 334, 312, 284 MeV, respectively. (A more sophis- 
ticated analysis will be performed in Sect. |v|.) 

We show, in Figs. |l2|, |l3|, the temporal correlators 
G(t)/G(0) together with the associated effective masses 
m eS (t) at T = 312,390 MeV, respectively, for the 0++ 
glueball-like color-singlet mode. In each effective-mass 
plot, a plateau is observed in the region indicated by the 
vertical dotted lines. By simply applying the single-cosh 



fit within the fit-range determined from the plateau, the 
"pole-mass" of the ++ glueball is extracted above T c , 
which is listed in Table B. We observe that the "pole- 
mass" is of about 1200 MeV, and that it changes rather 
continuously across the critical temperature T c . 

We perform the similar analysis for the 2 ++ glueball- 
like mode. In Figs. 14 and [ll| the temporal correlators 
G(t)/G(0) are plotted together with the associated effec- 
tive mass plots at T = 312,390 MeV, respectively. In 
each effective mass plot, we observe a plateau in the re- 
gion, which is indicated by the vertical dotted lines. Sim- 
ilar to the ++ case, we simply apply the single-cosh fit 
within the fit-range determined from the plateau to ex- 
tract the "pole-mass" . The results are listed in Table ||. 
We observe that the "pole-mass" is about 1400 MeV, 
and that its change is again rather continuous across the 
critical temperature T c . 

In both of these ++ and 2 ++ cases, the size of the 
plateau tends to shrink at high temperatures above T c . 
We note that, at much higher temperatures, the plateau 
finally disappears, which, however, seems to be due to 
the failure of the ground-state saturation in such a lim- 
ited temporal size. At any rate, these phenomena suggest 
that the single-cosh ansatz becomes inappropriate for the 
fit-function of the temporal correlator G(t) above T c at 
such a high temperature. However, because of the exis- 
tence of the plateau in the effective-mass plot, we do not 
exclude so far the possibility that the color-singlet modes 
exist as meta-stable modes just above T c . According to 
the lattice data, the color-singlet ++ and 2 ++ modes 
have the pole-mass of about 1200 MeV and 1400 MeV, 
respectively, provided that it can be regarded as a bound 
state. However, even in this case, it should be kept in 
mind that such glueball-like modes can decay into two or 
more gluons in the quark-gluon-plasma phase. Since they 
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TABLE II: The pole-mass ma(T) of the lowest 2 ++ glueball at finite temperature T in SU(3) lattice QCD. The meaning of 
T, Nt, vhg{T), (7 max i x 2 /N~df, (tiftz), g(0), Nsmi, R, and N con { are the same as those in Table |. The best smearing on iV smr 
is determined with maximum overlap condition. 



T[MeV] 


N t 


m G [MeV] 






(ti,t 2 ) 


9(0) 


iV m 


o[fml 

r L J 


7V CO nf 


130 


72 


2167(35) 


0.97(1) 


1.01 


( 5,13) 


0.97(1) 


54 


0.50 


5500 


187 


50 


2123(39) 


0.95(2) 


0.57 


( 5,13) 


0.95(2) 


39 


0.43 


5700 


208 


45 


2054(58) 


0.90(3) 


1.07 


( 7,19) 


0.90(3) 


47 


0.47 


6400 


218 


43 


2146(55) 


0.91(3) 


1.49 


( 7,12) 


0.91(3) 


68 


0.56 


9200 


234 


40 


1986(72) 


0.85(4) 


1.03 


( 8,14) 


0.85(4) 


63 


0.54 


8600 


246 


38 


2142(55) 


0.95(3) 


0.33 


( 6,14) 


0.95(3) 


41 


0.44 


8900 


253 


37 


2009(47) 


0.90(2) 


1.32 


( 6,13) 


0.90(2) 


38 


0.42 


8900 


260 


36 


1928(46) 


0.88(2) 


3.12 


( 6,12) 


0.88(2) 


52 


0.49 


9900 


268 


35 


1975(65) 


0.91(4) 


2.12 


( 8,16) 


0.91(4) 


44 


0.45 


9900 


275 


34 


1607(61) 


0.69(4) 


0.80 


( 9,14) 


0.69(4) 
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0.47 


9900 


284 


33 


1466(78) 


0.60(5) 


1.51 


(11,16) 


0.60(5) 
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1319(58) 


0.61(4) 


0.33 


(10,15) 


0.62(4) 
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1408(52) 
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1.93 


( 8,14) 
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0.47 
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1440(38) 


0.75(2) 


0.59 


( 7,13) 


0.76(2) 
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0.53 


6800 


390 


24 


1466(55) 


0.68(4) 


1.47 


( 8,12) 


0.70(3) 


32 


0.39 


7700 



cannot be stable above T c , they will acquire finite decay 
widths. Hence, for more reliable analysis, it is necessary 
to take into account the effect of such a width, which will 
be attempted in Sect. |y|. 

D. Discussions and remarks 

In this subsection, we compare our results of the pole- 
mass reduction of the thermal glueball with the related 
lattice QCD results on the screening masses as well as 
the pole-masses of various hadrons at finite temperature. 

We begin with the discussion of the pole-masses of 
mesons. The pole-masses of various mesons have been 
measured at finite temperature by using the anisotropic 
SU(3) lattice QCD at the quenched level in Refs. |J, (jj. 
From the same pole-mass analysis as ours, they have 
found that, both for light mesons and for heavy mesons 
such as the charmonium, the thermal pole-masses are 
almost unchanged from their zero-temperature values 
within the error bar below T c . This tendency persists 
even in the vicinity of T c as 



,(T = 0), for T<T C 



(44) 



In contrast, with the same pole-mass analysis, the lowest 
++ glueball shows the significant pole-mass reduction 
of about 300MeV in the vicinity of T c . Hence, the pole- 
mass reduction of the thermal ++ glueball may serve 
as a better pre-critical phenomenon of the QCD phase 
transition at finite temperature than the pole-mass shifts 
of the thermal mesons. Thus the thermal glueball may 
become an interesting particle in the RHIC project in the 
future. 

We next discuss the related results on the screening 
mass. The screening mass is defined as the reciprocal 
of the correlation length along the spatial direction. It 
can also reflect some of the important nonperturbative 



features of the QCD vacuum. Hence, it is also an inter- 
esting quantity, which may provide a signal of the QCD 
phase transition. However, unlike the pole-mass, which 
is considered to be directly measurable as a mass of an el- 
ementary excitation at finite temperature, the screening 
mass is not thought of as a directly measurable quantity 
in high energy experiments. In the ++ glueball sector, 
a significant reduction of the screening mass is reported 
as 



,: (T = 0.75T c )/mg r (r - 0) = 0.6 ± 0.1 



(45) 



in SU(3) isotropic lattice QCD with (3\ at — 5.93 over 
the lattice of the size 16 3 x N t with Nt — 4, 6, 8 at the 
quenched level |l^] . This significant reduction of the 
screening mass of the glueball would be an indication of 
the decrease of the nonperturbative nature of the QCD 
vacuum. In this sense, it seems to be consistent with our 
results on the significant reduction of the pole-mass of 
the thermal ++ glueball near the critical temperature 
T c . In contrast to the ++ glueball, in the meson sector, 
the changes of the screening mass has been found again 
rather small below T c at the quenched level Jf7|. 



V. NEW ANALYSIS OF TEMPORAL 
CORRELATIONS —SPECTRAL FUNCTION 

In this section, we attempt to generalize the pole-mass 
analysis given in Sect. 



[V. Our aim is to analyze our 



lattice QCD data of the temporal correlator G(t)/G(0) 
without assuming the narrowness of the thermal width 
of the bound-state peak in the spectral function. We 
first consider the functional form of the new fit-function, 
which takes into account the effect of the thermal width. 
We also provide its theoretical background. We then 
present the numerical results of this new analysis. The 
section is closed with several comments and cautions. 
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FIG. 9: (a) The 2 ++ glueball correlator G(t)/G(0) of the at 
the low temperature (T — 130 MeV). (b) The corresponding 
cosh- type effective-mass plot. The meaning of the solid lines 
in both figures and the dashed curve in (a) are the same as in 
Fig. ^ In (b), the points beyond t = 25 are suppressed due 
to their huge statistical errors. 



A. Theoretical consideration on the temporal 
correlator and the Breit-Wigner fit-function 



In Sect. |V|, we have presented the results of the pole- 
masses of the thermal ++ and 2 ++ glueballs at finite 
temperature obtained in the best-fit analysis of the tem- 
poral correlator G(t)/G(0) with the fit function of the 
single-cosh type. In this analysis, we have regarded 
each thermal glueball as a quasi-particle and have as- 
sumed that the thermal width of the corresponding peak 
is sufficiently narrow in the spectral function. In this 
case, the ground-state contribution to the spectral func- 
tion can be well approximated as p(u>) ~ 27r{<5(o; — 
mcC? 1 ))— S(uj+mQ(T))} by introducing the temperature- 
dependent pole-mass i7ig(T) of the thermal glueball. As 
a consequence, in order to measure the pole-mass, we can 
adopt the standard single hyperbolic cosine in Eq. J37| ) 
as the fit function, and perform the best-fit analysis in 



o 

| 




(b) 



2500 



2000 



> 

£ 1500 



1000 



„ 

' • • . 

















10 



15 



t [a,] 



FIG. 10: (a) The 2++ glueball correlator G(t)/G(0) at the 
high temperature (T = 253 MeV). (b) The corresponding 
cosh- type effective-mass plot. The meaning of the solid lines 
in both figures and the dashed curve in (a) are the same as 
in Fig. & 



exactly the same way as the zero-temperature case. 

However, to be strict, at finite temperature, each 
bound-state peak in the spectral function p(u>) acquires a 
non-zero width through the thermal fluctuations. Since 
the thermal width would grow with the temperature, it 
is expected to play the more significant roles and become 
the less negligible at the higher temperature. Hence, it 
is desirable to take into account the effect of the thermal 
width in the best-fit analysis. We can directly see how 
such a thermal width T(T) affects the temporal correlator 
G{t) in the spectral representation Eq. (g). 

There may have already appeared a signal of such 
a non-zero thermal width in the glueball correlators in 
Figs. ||, H, [l2| [l^, where the ground-state enhancement 
works the less effective for the higher temperature. In 
Table H, we observe a tendency that g(0) declines with 
the increasing temperature. For instance, g(0) ~ 0.95 for 
the low temperature T = 130 MeV becomes g(0) < 0.90 
for the higher temperature T > 250 MeV for the ++ 
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FIG. 11: The 2 ++ glueball mass plotted against temperature 
T. The solid circle and the cross denote the thermal 2 ++ 
glueball mass obtained with the maximum overlap condition 
and the mass minimum condition, respectively. The vertical 
dotted line indicates the critical temperature T c ~ 280 MeV 
in quenched QCD. 



glueball. In the effective-mass plots in Figs. 0, |[ |l^,we 
observe that the plateau tends to start with the larger t 
at the higher temperature. These phenomena may indi- 
cate that the existence of the other spectral components 
become the less negligible at the higher temperature. 

In fact, at finite temperature, each bound state peak 
acquires a thermal width. The thermal width is one of 
the most feasible candidates for the source of these other 
spectral components, since the contributions from such a 
bunch of spectral components cannot be well separated 
even by means of the smearing procedure. In this section, 
we perform a more general new analysis of the temporal 
correlator G(t)/G(0) as an attempt to take into account 
the effects of the thermal width T(T) of the ground-state 
peak in the spectral function p(uj). 

At zero temperature, the bound-state poles of the 
Green functions Gr(w) and G\{uj) in Eq. (||) lie on the 
real axis in the complex w-plane. Hence, the lowest-state 
contribution to the spectral function p{uj) is expressed as 
p(uj) — In A {5(u> — mo) — 5{uj + mc)}, where A and vtlq 
are the strength and the mass of the lowest-state pole, re- 
spectively. By inserting this expression into Eq. (0), and 
taking the limit (3 = 1/T — > oo, we recover the familiar 
expression as 



G(t) = lim A 

/3— »oc 



cosh (m G (/3/2 - t)) 
sinh(/3mG/2) 



= At 



—mat 



(46) 



With the increasing temperature, these bound-state 
poles are moving off the real axis into the second Rie- 
mannian sheet of the complex cj-plane. Suppose that the 
ground-state pole is located at uj = uj — £T (ui n , T e K) 
in the second Riemannian sheet. Since the spectral func- 
tion p(uj) is the imaginary part of the Green function, the 
ground-state contribution in the spectral function p(lo) 
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FIG. 12: (a) The ++ glueball correlator G{t)/G(0) at T = 
312 MeV above T c . (b) The corresponding cosh-type effective- 
mass plot. The meaning of the solid lines in both figures and 
the dashed curve in (a) are the same as in Fig. ^| 



can be expressed in the following way: 

pH = -2Im(G R ( W )) (47) 
= 27ivl (5r(u — wo) — S r {uj + o7 )) -I . 

Here, A represents the residue of the pole, S e (x) denotes 
the Lorentzian at x — with the width e (> 0), which is 
defined as 



8Jx) = — Im 

7T 



1 



1 



7T x 2 + e 2 



(48) 



The appearance of the second term in Eq. (|47|) is due to 
the fact that the spectral function p(u) is an odd func- 
tion in uj, reflecting the bosonic nature of the glueball. 
"• • •" in Eq. ( |47j ) represents the contributions from the 
excited states, which are expected to be suppressed after 
the appropriate smearing procedure. Corresponding to 
Eq. (fFfl); the ground state contribution in the temporal 
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FIG. 13: (a) The ++ glueball correlator G(i)/G(0) at T = 
390 MeV above T c . (b) The corresponding cosh- type effective- 
mass plot. The meaning of the solid lines in both figures and 
the dashed curve in (a) are the same as in Fig. |E| 



glueball correlator G(t) is expressed as 
du cosh(u;(/3/2 - *)) 



G(t) = 



,2tt sinh(/3u;/2) 

x2ttA (5r(u — wo) — <5r(^ + Wo)) + 



(49) 



Hence, to extract the center ujq{T) and the thermal width 
r(T) of the lowest-state peak in G(t)/G(0), the appro- 
priate fit-function is given as 



flr(*) 



do; cosh(a;(/?/2 - f)) 
3 2^ sinh(/3w/2) 
x2nA ^(5r(^ — ^o) — <5r(w + ^o)) 



(50) 



where A, o>o and the T are used as the fit parameters. 
Here, A corresponds to A/G(0) in Eq. (f49|), and will be 
referred to as the strength parameter. We will refer to 



Eq. (|5(j) as the fit-function of "Breit-Wigner" type. Note 

that gr(t) is a generalization of the fit-function g(t) of the 
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FIG. 14: (a) The 2 ++ glueball correlator G(t)/G(0) at T = 
312 MeV above T c . (b) The corresponding cosh-type effective- 
mass plot. The meaning of the solid lines in both figures and 
the dashed curve in (a) are the same as in Fig. ^| 



single-cosh type in the sense that gr(t) reduces to g(t) in 
a special limit as 



^^ffrti) = g{t). 



(51) 



In this sense, the best-fit analysis of Breit-Wigner type 
serves as a generalization of the ordinary pole-mass anal- 
ysis of single-cosh type. Note that the analysis of Breit- 
Wigner type is rather general, which is also applicable 
to the analysis of temporal correlators of various thermal 
hadrons. 

We note that, after some calculations, the fit-function 
gr(t) in Eq. ( pfj| ) can be expressed by the infinite series 
as 
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FIG. 15: (a) The 2 ++ glueball correlator G(t)/G(0) at T = 
390 MeV above T c . (b) The corresponding cosh- type effective- 
mass plot. The meaning of the solid lines in both figures and 
the dashed curve in (a) are the same as in Fig. |E| 



gr(t) = A 



Re 



cosh{(o;o + «r)£} 



sinh 



2/3^0 



OO 

E 

n=l 



COS 



27m 



1 



(2im + f3T) 2 + [3 2 



-(n 



(52) 



We make a comment on the role of the smearing 
method in the presence of the non-zero width of the low- 
est state (or the ground state) at finite temperature. As 
was explained in Sect. [Eij, the smearing method aims at 
a suitable choice of the glueball operators, by providing 
a series of operators with different sizes, all of which, 
however, hold the identical quantum numbers in com- 
mon. Whereas the residues of the bound-state poles are 
affected by the different choices of the glueball operators, 
their positions in the complex w-plane are not affected. 
One can thereby enhance the residue of the particular 
bound-state pole by a clever choice of the operator, while 



keeping its position unaffected. Note that the position of 
the pole determines the characteristics of the peak, i.e., 
the center and the width. As this clever choice of the 
operator, we will use the smearing method, which is ex- 
pected to pick up and then to enhance the contribution 
of the ground-state peak. 



B. Setup for the Breit-Wigner fit analysis 

In this subsection, we consider the fit-range of the 
Breit-Wigner fit analysis for the temporal glueball cor- 
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(b) T = 253 MeV 
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FIG. 16: The effective center ojo;off(i) and the effective width 
r e ff (t) of the ++ glueball correlator at various temperatures 
(a) T = 130 MeV, (b) T = 253 MeV and (c) T = 390 MeV. 
The solid lines represents the results of the best-fit analysis 
of Breit-Wigner type. Here, the fit-ranges are determined 
from the simultaneous plateaus, which are indicated with the 
vertical dotted lines. 



relator G(t). We examine also the iV smr -dependence of 
the best-fit parameters in order to estimate the system- 
atic error originating from a particular choice of N BmT . 



1. The fit range 

We consider the determination of the fit-range of the 
Breit-Wigner fit analysis. As a consequence of the smear- 
ing method, the contribution from the ground-state peak 
is expected to be enhanced in the spectral function p(to). 
However, even with the smearing method, the complete 
elimination of all the excited-state contributions is prac- 
tically impossible. Therefore, we have to reduce fur- 
ther the remaining contributions from the higher spec- 
tral components as much as possible. To this end, we 
seek for the appropriate fit-range, where G(t) consists of 
nearly a single-peak contribution. We adopt the strategy, 
which is a straightforward extension to the one adopted 
in the single-cosh analysis, i.e., the analysis based on the 
effective-mass plot in the following way. We consider 
the "effective center" a>o; C ff(i) and the "effective width" 
r c ff(i), which are defined as the solutions to the following 
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FIG. 17: The normalized overlap gr(0), the center uio and 
the thermal width F of the lowest peak of the ++ glueball 
against the smearing number N aulI in the confinement phase 
at T = 130 MeV and 250 MeV, which are denoted by the circle 
and the triangle, respectively. The solid symbols denote the 
data where <?r(0) becomes maximum. Both for u>o and T, the 
iV smr -dependence is small in the region of 30 < N SII1I < 50. 



equations: 



G(t)/G(t + 1) 
G(t + l)/G(t + 2) 



g r (t)/gr(t + l) (53) 
gr(t + l)/gr(t + 2), 



for given (M/G(t+1) and G(t+l)/G(t+2) at each fixed 
t. In Fig. Ilfl, we show the plots of the effective centers 
i^o-,es{t) and the effective width r e gf(i) of the ++ glue- 
ball correlators at various temperatures T = 130, 253, 390 
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FIG. 18: The similar figure as Fig. [l?] in the deconfinement 
phase at T = 312 MeV and 390 MeV, denoted by the square 
and diamond, respectively. The absence of the solid diamond 
is due to the missing maximum of the normalized overlap 
gr{0) in the region N Bml < 70. 
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FIG. 19: The spectral function p(ui) of the lowest ++ glue- 
ball at T = 130, 253, 312, 390 MeV, obtained from the Breit- 
Wigner fit analysis of the ++ -glueball temporal correlator 
G(t)/G(t). 
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FIG. 20: The spectral function p(ui) of the lowest 2 ++ glue- 
ball at T = 130, 253, 312, 390 MeV, obtained from the Breit- 
Wigner fit analysis of the 2 ++ -glueball temporal correlator 
G(t)/G(t). 



of the simultaneous plateau in the pair of coo- e g(t) and 
r c ff(i) plots works as a necessary condition to determine 
whether the fit function gr(t) of Breit-Wigner type is 
appropriate or not. 



MeV. The statistical errors for wo ;e ff(i) and T c g(t) are 
estimated with the jackknife analysis. In each pair of 
wo;eff(i) an d F e ff(£) plots, there appears a simultaneous 
plateau, a region where ujQ- c s(t) and T e g(t) are almost 
constant simultaneously. In this region, G(t) is expected 
to consist of a single peak contribution, and grit) can 
properly represent the contribution of the ground-state 
peak in Eq. (49). Note that a nontrivial t-dependence 



indicates the existence of the contributions from the 
other spectral components. In this sense, the existence 



2. The N s , 



We consider the A^ smr -dependences of the best-fit pa- 
rameters, the center ojo(T) and the thermal width T(T) 
of the lowest peak. Since the physical meaning of A is 
not so obvious in the presence of the thermal width, we 
consider gr(0) instead of A. Note that <?r(0) is a gener- 
alization of g(0) in the single-cosh analysis. We will refer 
to <?r(0) as the "(normalized) overlap" with the state cor- 
responding to the lowest peak. 
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TABLE III: The analysis of the temporal correlator G(t)/G(0) of the thermal ++ glueball based on Breit-Wigner ansatz of 
the spectral function p(uj) as Eq. mt). The temperature T, the temporal lattice size Nt, the center coo, the width F and the 
strength parameter A of the peak, the correlated x /Ndf, the fit-range (ti,ta)) the overlap <7r(0), the smearing number iV smr 
which achieves the gr(0)-max condition, and the gauge configuration number iV con f are listed. The asterisk "*" in the column 
of JVsmr indicates that the original value of JV smr is 70, which had been replaced by the suitable smearing number 40. We had 
done this, because N eml = 70 is the end point and because the situation similar to the one in Fig. hq at T = 390 is observed. 
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TABLE IV: The analysis of the thermal 2 ++ glueball correlator G(t)/G(0) with Breit-Wigner ansatz of the spectral function 
p(uj). The meaning of T, Nt, ujq, V, A, x 2 /Ndf, (ti,ta)) ffr(0) and A^conf are the same as those in Table Eli . 
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In Fig. [l?], the normalized overlap <?r(0), the center 
u>o(T) and the thermal width T(T) are plotted against 
N smI for fixed a = 2.1 at T = 130 and 253 MeV below 
T c with circles and triangles, respectively. In the region of 
30 < N sml < 50, flfr(O) is maximized, and lu (T) and T(T) 
are minimized. These behaviors of gr(0) and loq(T) are 
analogous to those of the normalized overlap g(0) (or the 
overlap C) and the pole- mass mc(T) in the single-cosh 
analysis. As for the thermal width T(T), since the con- 
tamination of the higher spectral components is expected 
to make it wider, the minimal T(T) would characterize 
the lowest-peak saturation. 

According to these considerations, we expect that the 
lowest-peak contribution is maximally enhanced in the 



region of 30 < N smI < 50, and that the normalized over- 
lap gr(0) — 1 provides a saturation rate of the lowest- 
peak contribution, in exactly the same way as g(0) ~ 1 
does in the single-cosh analysis. Note that both for the 
center coq and the thermal width T, the A^ sinr -dependence 
is rather small in the region of 30 < N sml < 50, result- 
ing in the small systematic error on the specific choice 
of N smI . The behaviors are qualitatively the same at the 
other temperatures below T c . This is also the case for 
the 2++ glueball. 

In contrast, qualitatively different behaviors appear 
above T c . In Fig. [l8|, the normalized overlap <7r(0), the 
center u>o and the thermal width T are plotted against 
N smI for a = 2.1 at T = 312 and 390 MeV above T c 
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FIG. 21: The center oj (T) and the width Y{T) of the lowest- 
state peak in the spectral function p{ui) for the ++ glue- 
ball against temperature T in the Breit-Wigner fit analysis. 
The triangle denotes the results associated with the suitable 
IK and the circle denote the results with 
Ll]L which maximizes the normalized over- 



smearing as Eq. (2' 
iV amr listed in Tab 
lap g r (0). 
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FIG. 22: The center w (T) and the width T(T) of the lowest- 
state peak in the spectral function p{u) for the 2 ++ glue- 
ball against temperature T in the Breit-Wigner fit analysis. 
The triangle denotes the results in the suitable smearing as 
Eq. (I29T), and the circle denote the results with N smT listed in 
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Table IV, which maximizes the normalized overlap <?r(0). 



3. The smearing and the single-pole saturation 



with the square and the diamond, respectively. We see 
that the overlap gr(0) becomes too insensitive to iV smr to 
determine its maximum. In some cases, it does not take 
the maximum in the region iV smr < 70. In addition, there 
are sizable -/V smr -dependences in loq and T in this region, 
although r seem to take its minimum around the typical 
smearing number as -/V smr ~ 40. As for the center loq, 
it continues to decline and does not take the minimum 
in the region iV smr < 70. We note that qualitatively the 
same behaviors are observed in the 2 ++ channel. Consid- 
ering these behaviors, it would be better to consider the 
Breit-Wigner analysis both in the case with N amT = 40 
and in the case with N smT which maximizes gr(0) above 



In this section, we use the smearing method to enhance 
the residues of the pair of the complex poles associated 
with the low-lying peak. However, in the presence of the 
non-zero thermal width, the meaning of the "low-lying" 
contribution may become less definite in a strict sense, 
since these peaks can overlap with each other. With 
the increasing temperature, the thermal width becomes 
wider. Hence, each bound-state peak becomes less dis- 
tinguishable from one another due to the possibly strong 
overlap with the neighboring peaks. In some cases, one 
may wonder if the smearing really extract only the con- 
tribution from the lowest-state peak among the overlap- 
ping peaks. Therefore, we have to check the smallness of 
the excited-state contamination by examining the qual- 
ity of the fitting of G(t)/G(0) with g T {t). If G(t)/G(0) 
can be well fitted with gr(t) corresponding to the sin- 
gle Lorentzian peak, we may expect that the lowest-state 



23 



dominates in the spectral function p(ui) as a result of the 
smearing. 

C. Numerical results of the Breit-Wigner fit 

In this subsection, we perform the Breit-Wigner fit 
analysis for the lattice QCD data of the temporal glue- 
ball correlator. We also present the numerical results on 
the spectral function of the glueballs at various temper- 
atures. 

1. The glueball correlators and the spectral functions 



In Figs. f[ 



), P (a), 


12 (a) 




13 (a) 



we show the 



best- fit function gr{t) for the ++ glueball correlator 
G(t)/G(0) with dashed curves at various temperatures 
as T = 130,253,312,390 M eV, r e spectively. S i milarly 
for the 2++ glueball, in Figs. |9 (a)|, [LP (a% |14 (a)|, [15 (aj 



the best- fit functions gr (t) are added with dashed curves 
at various temperatures as T = 130,253,312,390 MeV, 
respectively. It is remarkable that, at all of these tem- 
peratures, the lattice QCD data of the temporal glueball 
correlators G(t) /G(0) are in a good agreement with gr(t) 
in the whole region of t. 

In Figs. |l^ and [2(], we plot the spectral functions p(uj) 
in the ++ and 2 ++ glueball channels at various temper- 
atures. In both the channels, we observe the tendency 
that the thermal width grows with the increasing tem- 
perature (up to the statistical errors). To be strict, each 
of these curves corresponds to actually a part of the spec- 
tral function p(uj), i.e., it is the ground-state contribution 
in the spectral function rather than the original spectral 
function itself. This can be understood in the follow- 
ing way. Recall that the best-fit analysis is performed in 
the fit-range determined from the simultaneous plateau 
in the plots of the effective center wo;eff(i) and the effec- 
tive width r c ff(£). In this fit-range, G(t) is expected to 
be less sensitive to the existence of the higher spectral 
components, and hence, it becomes possible to extract 
only the contribution from the ground-state peak. How- 
ever, in order to express G(t) with the spectral function 
p{uj) (cf. Eq. (||)) in the whole region of t, the higher 
spectral components should be necessary even after the 
smearing procedure is applied. Nevertheless, the good- 
ness of the fitting in Figs. |, |, [jj, [jj, §, [h], || and |l| 
indicates that the contributions from the higher spectral 
components may be actually rather small. 



2. The main result of the Breit- Wigner fit 

In Fig. |2l], the center ujq(T) and the thermal width 
T(T) of the lowest-state peak of the spectral function 
p(u) in the ++ glueball channel are plotted against tem- 
perature T. Here, the triangles denote the results for the 
suitable smearing with N sml = 40, and the circles denote 



the results with N smr , which maximizes the normalized 
overlap #r(0). The agreement of both results indicates 
that the systematic error originating from the particular 
choice of N smx is small. 

We observe that the thermal width T(T) grows grad- 
ually with the increasing temperature, resulting in the 
significant increase near the critical temperature T c as 



T(T ~ T c ) ~ 300MeV. 



(54) 



The center luq(T) is observed to decline modestly by 
about 100 MeV near T c as u Q {T ~ T c ) ~ 1440 MeV 
in comparison with u>o(T ~ 0) ~ 1530 MeV. 

In Fig. [2^, the center ujq(T) and the thermal width 
T(T) of the lowest-state peak of the spectral function 
p(uj) in the 2 ++ glueball channel are plotted against tem- 
perature T. In comparison with the ++ glueball case, 
changes are observed less significant in the 2 ++ glueball 
channel below T c , which, however, is followed by a sud- 
den expansion of the thermal width T(T) around T ~ T c . 
This sudden expansion of T(T) may indicate the insta- 
bility of the thermal 2 ++ glueball around T c . 

In Tables [nj| and [Fv|, we summarize the results of the 
Breit-Wigner fit analysis for the ++ and the 2 ++ glue- 
balls, respectively, such as the center uq(T), the thermal 
width r(T) and the strength parameter A of the lowest- 
state peak, which play the role of the fit parameters, and 
so on. 

As a remarkable result of the Breit-Wigner analysis, 
we emphasize that the gradual growth of the thermal 
width r(T) of the ++ glueball begins already far below 
T c in the confinement phase. This may be an attractive 
feature for the experimental observation of the thermal 
effect. 



3. Comparison of the Breit- Wigner analysis with the 
pole-mass analysis 

Finally, we consider the relations between the pole- 
mass analysis and the Breit-Wigner analysis. As is men- 
tioned before, the Breit-Wigner analysis is a straightfor- 
ward extension to the pole-mass analysis. In fact, they 
coincide in the vanishing width limit r — * +0. In this 
limit, we have the equality 



m G {T) = u) (T). 



(55) 



However, in the presence of the width T, the relation 
between ma(T) and co$(T) becomes nontrivial, and we 
are left with the inequality, 



m G {T) < u {T). 



(56) 



To see this, we consider the spectral representation 
Eq. (||), where the temporal correlator G(t) is expressed 
as an average of the hyperbolic cosine of frequency u> with 
the weight as 



2sinh(/3w/2)' 



(57) 
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Suppose the idealized case where the spectral function 
p(tjj) is completely dominated by a single peak. In this 
case, Eq. ( |57| ) reads 

2-rtA [Sr(u> — uq) — Sr(oJ + ojq) ) 
W M= 2sinh W 2) ^ 

To define the pole-mass, we approximate the weight 
W(uj) with the (5-functions as 



27ivl (s(uj - m G ) + S(lo + m G )j 



2sinh(/3m G /2) 



(59) 



where the pole-mass tog is determined so as to reproduce 
the shape of G(i)/G(0) in a fit range as faithful as possi- 
ble. Note that sinh(/3a->/2) in the denominator of Eq. d5g| ) 
works as a biased factor, which enhances the smaller lo 
region while suppressing the larger u> region. As a conse- 
quence of this biased factor, the peak position of W(ui) 
is shifted from ujq to a smaller value, which leads to the 
inequality Eq. (|5^). 

To make a rough estimate of the pole-mass tog, we 
consider the peak position uj = cj max (> 0) of the weight 
W(lu) in Eq. (|58|), which is defined as 



du> 



W(uj) 



= 0. 



We adopt the simple identification as 
fiG — w max . 



(60) 



(61) 



In the confinement phase T < T c , the results of the Breit- 
Wigner fit indicate the inequality as 

w >r,T. (62) 

In this case, Eq. ([58]) can be approximated as 

W(lu) ~2TrAe~ ,3u ' /2 5r(uJ ~lu ) (63) 

in the region of w ~ ujq. Hence, Eq. (pG) reduces to 



uj ) 2 + r 2 } ~o. 



(64) 



As a consequence of the simple identification in Eq. (61 ), 
we find the relation among toq(T), u)q(T) and T(T) as 

to g (T) ~ wo(T) - {2T - V4T 2 -r(T) 2 } < wo(T), 

(65) 

which seems consistent with our numerical results below 
T c . 

We make a comment on the experimentally observed 
particle mass at finite temperature in the presence of the 
width r. In the high-energy experiments, the spectral 
function p(u>) is actually observed, providing the relevant 
information of the mass and the width of the particle at 
finite temperature. In fact, the experimentally observed 



mass is expected to be distributed around to = u)q(T) 
with the width T(T). As long as the width T is enough 
narrow, the pole- mass tog(T) provides a good approx- 
imation of the peak position of l<Jq(T) in p(co). In this 
case, the pole- mass itlq(T) can be regarded as the ther- 
mal particle mass. In fact, in all of the previous stud- 
ies of the hadronic temporal correlations |l^, [l4|, fl6| , 
the single-cosh analysis has been used assuming that the 
thermal width is enough narrow. However, when the 
width r becomes wide, the pole-mass toq(T) obtained 
from the single-cosh analysis suffers from the effect of 
the biased factor "sinh(/3o;/2)" in the denominator in 
Eq. (|57|). Therefore, tog(T) is shifted from ujq to a 
smaller value. In this case, the Breit-Wigner analysis is 
preferable. As is demonstrated before, the Breit-Wigner 
analysis physically means the direct measurement of the 
pole position of the low-lying particle in the complex em- 
plane, i.e., the shape of the low- lying peak in the spectral 
function p(u>). It provides the mass and the width of the 
particle at finite temperature as the center u)q(T) and the 
width r(T), respectively. 



VI. SUMMARY AND CONCLUDING 
REMARKS 

We have studied the thermal properties of the ++ 
and 2 ++ glueballs using SU(3) anisotropic lattice QCD 
with /3i at = 6.25, the renormalized anisotropy £ = 
a s /a t = 4 (a„ ~ 0.084 fm and a t ~ 0.021 fm) 
over the lattice of the size 20 3 x N t with N t — 
24,26,28,30,33,34,35,36,37,38,40,43,45,50,72 at the 
quenched level. 

To begin with, we have measured the temporal corre- 
lators G(t) for the lowest ++ and 2 ++ glueballs using 
more than 5,500 gauge configurations at each tempera- 
ture. In this calculation, we have adopted the smearing 
method to construct the suitable operator for the lowest- 
lying ++ and 2 ++ glueballs on the lattice. We have 
also provided an analytical consideration on the physical 
meaning of the smearing procedure based on the spatial 
distribution of the gluon field in the smeared operator. 

Next, we have performed the pole-mass measurements 
of the thermal glueballs from G(t) by adopting the pro- 
cedure used in the standard hadron mass measurements. 
For the lowest ++ glueball, we have observed a signif- 
icant pole-mass reduction of about 300 MeV near T c or 
tog(T ~ T c ) ~ 0.8tog(T ~ 0), while its size remains 
almost unchanged as p(T) ~ 0.4 fm. This pole-mass 
shift is actually much larger than any other pole-mass 
shifts which have been ever observed in the meson sector 
fl3| |l4| with the similar analysis in lattice QCD. 

Finally, for completeness, we have performed a more 
general new fit-analysis of the temporal glueball correla- 
tor G(t) as an attempt to take into account the effects 
of the possible appearance of the thermal width T(T) 
of the bound-state peak. We have proposed the Breit- 
Wigner form as a generalized fit- function gr(t) for the 
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lowest-peak in the spectral function p(u>) of the temporal 
correlator G(t) at finite temperature. This ansatz is also 
applicable to the temporal correlators of various ther- 
mal hadrons. In this advanced analysis, we have found 
a significant broadening of the lowest-glueball peak as 
r(T) ~ 300 MeV near T c as well as a rather modest 
reduction in its center of about 100 MeV. 

We have investigated also the temporal correlators of 
the color-singlet modes corresponding to the glueballs in 
the deconfinement phase. We have found that the ther- 
mal width r(T) increases monotonically with increasing 
temperature both in the ++ and the 2 ++ glueball chan- 
nels. However, from only our lattice data, it is difficult 
to determine whether such color singlet modes really sur- 
vive above T c or not as a meta stable modes, which is left 
for the future studies. 

In this way, these two analyses have indicated the sig- 
nificant thermal effects to the lowlying thermal ++ glue- 
ball near T c , such as the considerable pole-mass reduc- 
tion or the width broadening. Note that the actual decay 
width of the glueball candidates, /q(1500) and /q(1710), 



are known to be about 100 MeV, which is rather small. 
Hence, either the 300 MeV pole-mass reduction or the 
thermal width broadening of 300 MeV, if happens, will 
be seen in the change of the spectral function in the high- 
energy experiment. Therefore, the thermal properties of 
glueball may provide an interesting signal of the precrit- 
ical phenomenon of the QGP creation in the future ex- 
periment in RHIC. For more direct comparison with the 
experiment, it is desirable to include the dynamical-quark 
effect. 
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